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Abstract 



Previously developed Pade-related method of resummation for QCD ob- 
servables, which achieves exact renormalization-scale-invariance, is extended 
so that the scheme-invariance is obtained as well. The dependence on the 
leading scheme parameter C2 is eliminated by a variant of the method of the 
principle of minimal sensitivity. The subleading parameter C3 in the approx- 
imant is then fixed in such a way that the correct known location of the 
leading infrared renormalon pole is reproduced. Thus, /3-functions which go 
beyond the last perturbatively calculated order in the observable are used. 
The /3-functions in the approximant are quasianalytically continued by Pade 
approximants. Two aspects of nonperturbative physics are accounted for in 
the presented resummation: a mechanism of quasianalytic continuation from 
the weak- into the strong-coupling regime, and the (approximant-specific) 
contribution of the leading infrared renormalon. The case of the Bjorken po- 
larized sum rule is considered as a specific example of how the method works. 
PACS number(s): ll.10.Hi, 11.80.Fv, 12.38.Bx, 12.38. Cy 



I. INTRODUCTION 



In QCD, as a result of extensive perturbative calculations, some observables are now 
known to the next-to-next-to-leading order (NNLO, ~a 3 ) in the power expansion in the 
strong coupling parameter a = a s /Tr. Knowing such truncated perturbation series (TPS), 
the question of their resummation is gaining importance, especially if the typical process 
energies associated with the observable are low and thus the relevant coupling parameter is 
large. In such cases, it is to be expected that additional perturbative and nonperturbative 
effects, not explicitly contained in the TPS, will be numerically important. Many methods 
of resummation, based on the available TPS, try to incorporate such effects. Some of these 
methods eliminate the dependence on the renormalization scale (RScl) and scheme (RSch) 
by fixing them in the TPS itself in a judicious way - these methods could be regarded as 
renormalization-group-improved methods of resummation: BLM-fixing motivated by the 
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large-n/ considerations |T[], Stevenson's principle of minimal sensitivity (PMS) 0, Grun- 
berg's effective charge method (ECH) || (cf. Ref. [|J for a related method). Some of the 
recent works on resummations are the method of "commensurate scale relations" and related 
approaches || , a method using an analytic form of the coupling parameter || , ECH-related 
approaches 0, and an expansion in the two-loop coupling parameter ||. In the past few 
years, Pade approximants (PA's) were also shown to be a rather successful method of re- 
summation ||, especially since the resummed results show in general weakened RScl- and 
RSch-dependence. The diagonal Pade approximants (dPA's) are particularly well motivated 
for observables since they are RScl-independent in the approximation of the one-loop evo- 
lution of the coupling a s (Q 2 ) |nj. In addition, PA's go in their form beyond the polynomial 
form of the TPS on which they are based, thus contain a strong mechanism of quasiana- 
lytic continuation from the weak- into the strong-coupling regime, and can consequently 
incorporate some of the nonperturbative effects into the resummed result. 

Recently, an extension of the method of dPA's was presented |TTJ which leads to the exact 



perturbative RScl-independence of the resummed result. We then extended the method so 
that it is applicable also to NNLO TPS's fl2|| , and suggested there a way to fix the RSch 
by applying the principle of minimal sensitivity (PMS). The way suggested in |12| does not 
work properly in practice since no minimum of the PMS equation dA/dc2 = [Eq. (40) 
there] exists. The dependence of our approximants on the RSch-parameters c 2 and C3 of the 
original TPS is a significant problem when the approximants are applied to the low-energy 
observables like the Bjorken polarized sum rule (BjPSR) at low Q P hoton « GeV [0. 

This problem is addressed in the present paper. For the case of NNLO TPS, an ex- 
tended version A of our approximant is constructed where the dependence on the leading 
RSch-parameter c 2 is eliminated by a variant of the PMS. Subsequently, the sub-leading 
RSch-parameter c 3 is adjusted so that the approximant reproduces the correct location of 
the leading infrared renormalon pole. The latter procedure is carried out in the concrete ex- 
ample of the BjPSR. The same method of C3-fixing is then applied to Grunberg's ECH and 
Stevenson's TPS-PMS approximants. Hence, in the approximants we use /3-functions which 
go beyond the last perturbatively calculated order in the observable (NNLO). Further, a PA- 
type of quasianalytic continuation for the /3-functions is used in all these approximants. The 
resulting predictions for a^ s (M|) from the BjPSR are presented, along with those when 
PA's are applied to the BjPSR, and compared with the world average. Differences between 
between our approximant and the other methods (ECH, TPS-PMS; PA's) are pointed out. 

II. CONSTRUCTION OF C 2 INDEPENDENT APPROXIMANTS 

Consider a QCD observable S with negligible mass effects and known NNLO TPS 

S[ 2 ] = a (l + nao + r 2 al) , (1) 
with: a = a(ln Ql; c { 2 } , cf } , ■ • •) , n = n (In Qq) , r 2 = r 2 (ln Q 2 ; c { 2 0) ) . (2) 

We denoted a = a s /ir; Qq is the Euclidean RScl; cf (j > 2) are the RSch-parameters used 
in the TPS. The coupling parameter a(\nQ 2 ; . . .) in this RSch evolves according to the 
renormalization group equation (RGE) da/d\n.Q 2 = (3(a; c^, . . .). Here, the /3-function 
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has the power expansion (3(a) = — (3oa 2 (l+cia+C2^a 2 +cf^a 3 + ■ ■ •), and f3 and C\ are RScl- 
and RSch-invariant. This RGE can be integrated (see Appendix of 0]) 



ft In [|f")=- + Cl ln 
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where a = dQnQ 2 ,; c 2 , cjj , . . .) and A is a universal scale (~0.1 GeV). When subtracting 
(||) from the analogous equation for a = d(lnQ 2 ; c 2 , c 3 , . . .), an equation is obtained which 
relates a with do, i.e., determines a in terms of do- This equation then determines also the 
expansion of a in powers of do- From now on, we fix the "A-convention" to A = A. 

We make the following ansatz for our approximant, motivated by the RScl-invariant (but 
not RSch-invariant) approximant of Refs. [I2"],IB]: 
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... (j = 1, 2) as fixed numbers, and Co 



where we regard now the parameters c 2 . , c 3 

Five parameters in the approximant (d, Q\, Q\, c^ 1 ', c^') can be fixed by applying five 
conditions to the approximant. Three conditions are obtained from the so called minimal 
requirement: When we expand the approximant back in powers of do, the first three coef- 
ficients of the original TPS (P have to be reproduced. The additional two conditions are 
obtained by a variant of the PMS 
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This allows us to fix ci , s. If we took in 



c (1) 

c 2 



-c (2) 

' c 2 



'2 



I , and Cfc 



.(2) 



(5) 

(A; > 3), i.e., the 

approximants of we would obtain dAs^jdc 2 = — lOcido + C^do) ^ a®, i.e., the PMS 

condition would not be satisfied. This is the main reason why we take two different (leading) 
parameters c 2 in the two d's in (ffl). Further, since the two energy scales in (f|), and in the 



approximants of (T^|T^], are Q\ ^Q\ y it does not appear unnatural to have t^c 2 ■ But 



the (subleading) parameters c 3 cannot be fixed by such an approach since 



(dA s s/dc^)\ Sc3 = 2al + O(a« ) where : cW = (#'+#>)/2, 6<* = (#>-£>) . (6) 
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The same problem arises in the NNLO polynomial approximants ECH and TPS-PMS where 
dA s /dc 3 = (l/2)a* + O(a 5 ) d(As) 2 /dc 3 = a{j + 0(ag)]. We will take, for simplicity, 
c 3 =c 3 2 ' ) = c 3 , and the value of c 3 will be fixed later. 

Conditions (H) then depend also on 604 = (C4 ^ C4 ) which we set equal to zero to 
avoid further (presumably unnecessary) complications. Then the set of the five equations 
determining d, Q 2 and (.7 = 1,2) reads 
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Here, p 2 is an RScl- and RSch-invariant, and therefore it is straightforward to see that 
the solutions of the system flTD-flTDP for Q 2 and c 2 (j = 1, 2) and for a are independent of 

the original choice of the RScl (Ql) and of the RSch (c 2 °^ , C3 , . . .). Eqs. originate 
from PMS conditions (|5p, and the other three identities from the minimal condition. In 
particular, the latter three identities (|7[) and (pi]) show that a and the scales Q\ and Q\ are 
(^-independent irrespective of whether 5c 2 ^0 or 5c 2 = 0. 

The coupled system of three equations (|7|)~(0) f° r the three unknowns c 2 (j = l,2) and 
y_ = p Q ln(Qi/Q 2 ) can be solved numerically. The solutions which give |a| <C 1 or \a\ ;§> 1 
must be discarded because they would cause numerical instabilities in the approximant, and 
they would not make sense physically either - one of the scales Qi, Q 2 would be orders of 
magnitude different from the other. There are apparently two possibilities: 1.) c 2 and 
5c 2 are all real numbers; 2.) c| is real, y_ and 5c 2 are imaginary numbers. In both cases, 
the approximant itself would be real, as it shoud be. 

If there are several solutions which give different values for the approximant, we should 
choose (again within the PMS-logic) among them the solution with the smallest curvature 
with respect to c 2 and c 2 2 \ 



III. APPLICATION TO THE BJORKEN POLARIZED SUM RULE; C3-FIXING 

The Bjorken polarized sum rule (BjPSR) involves the isotriplet combination of the first 
moments over xbj of proton and neutron polarized structure functions 
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where p 2 = —Qp^< is 7* momentum transfer. At Q 2 h = 3GeV 2 where three quarks are 
assumed active (ri/ = 3), and if taking MS RSch and RScl Q 2 ] = Q 2 h , we have |l4] , p!5| : 
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Solving numerically the system of equations ©-(0), we obtain one solution only 

4 1} = 1.465 , 4 2) = 5.137 , Qi = 0.594 GeV , Q 2 = 1.164 GeV (=>- a = 0.3301) . (16) 

This solution is independent of the choice of RScl and RSch. For the time being, we will 
set the higher parameters cj^ = (k > 4, j — 1, 2). Now our approximant depends only on 
the still free parameter C3. This dependence is numerically significant. For a typical value 
a = 0.09 [=>■ ap(3GeV 2 )^ 0.283, ap(Mf) «0.113], the approximant (§ gives 0.1523 and 
0.1632 when C3 = 0,c 3 vIS , respectively, i.e. a difference of 7.2%. In the case of the ECH and 
TPS-PMS approximants for the BjPSR, the respective differences are 3.8% and 4.0%. 

Parameter C3 characterizes the N 3 LO term in the corresponding /3-functions [cf. Eq. (Q)], 
and information on its value cannot be obtained from the NNLO TPS [cf. Eq. (|6])]. Therefore, 
to fix C3, we should incorporate into the approximants a known piece of (nonperturbative) 
information beyond the NNLO TPS ([14]) . Natural candidates for this are the known locations 
of the poles fli"6l|i"7| of the leading infrared renormalon (IRi: z po \ c = l//?o) or ultraviolet 
renormalon (UVi: z po \ c = —1 / j3 ) , i.e., the poles of the Borel transform Bg(z) of S closest to 
the origin. Large-/3 evaluations |17||, based on the formulas of [[L(J and using simple Borel 
transforms in a variant of the V-scheme (RScl Q = Q p h exp(— 5/6) and one-loop-evolved 
a), suggest that the UVi contributions to the BjPSR at Qp h = 2-3 GeV 2 are suppressed in 
comparison to the IRi contributions by a factor 3-4 (cf. their Fig. 2). 

Therefore, we will fix C3 in the three approximants by incorporating in them the infor- 
mation on the location of the IRi pole z po \ e = l/Po (=4/9). For that, we employ RScl- and 
RSch-invariant Borel transforms. Simple Borel transforms are not RScl/RSch-invariant, 
the use of their TPS's leads to RScl/RSch-dependent c 3 -fixing, which we want to avoid. 
We use a variant of the invariant Borel transform B(z) introduced by Grunberg [IB , who 



in turn introduced it on the basis of the modified Borel transform of Ref. [19] 



r°° r 1 

S(Q 2 ^) = I dzexp [- Pl (Q 2 ph )z\ B s (z) , (17) 



where Pl is the first RScl/RSch-invariant § of S: Pl = -n+p ln(Q 2 /A 2 ) = O ln(Q 2 h /A 2 ). 
Here, A is the universal scale of Eq. (|3]), and A a scale which depends on the choice of S but is 
RScl/RSch-invariant and Q p h~ m dependent. The Pi(Q p h) is, up to an additive constant (the 



latter not affecting the poles of Bs), equal to 1/a^ 1 loop ' (QL) . Thus, B s (z) of dl7f) reduces 



to the simple Borel transform, up to a factor exp(cz), once higher than one-loop effects are 
ignored. The coefficients of the power expansion of B$(z) of (|T7|) are RScl/RSch-invariant, 
in contrast to the case of the simple Borel transform. These invariant coefficients can be 
related with coefficients r n of S most easily in a specific RSch Ck = c\ (k > 2) 



V (fn — Cifn-l) 



B s (z) = (c lZ r* eM-m) E r(n + i + c^) zW = t ^ 13 ^ " ( 18 ) 

Here, r n is the coefficient at a n+1 in the expansion of S in powers of a = a(lnQQ; c 2 , cf, ...); 
by definition f_i = 0, fo = 1. Thus, the expansion of the approximant y/ As^{cz) in powers 
of a leads to the expansion of the (reduced) Borel transform B^(z) in powers of z. The 
coefficients starting at z 3 are predictions of the approximant and c 3 -dependent: Br^(z) = 
l+b lZ +b 2 z 2 + b 3 z 3 +---, with 6x^-0.7516, b 2 ^0.4209, 63 w (-2.664+0.1667c 3 ), etc. Terms 
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with high powers of z are not reliable, because the approximant is based on an NNLO TPS 
S[2] with only two terms beyond the leading order. We then employ Pade approximants 
(PA's) of power expansion of B^r^, since they are efficient in determining the pole structure 
of Br%. We performed the expansion of v^4(c 3 ) up to ~a 7 , obtaining the expansion of 
Brj(z) up to ~^ 6 . This allowed us to construct PA-g's of as high order as [3/3] or [4/2]. 
The value of c 3 in PA-g was then adjusted to achieve z po i e = l//3 (=4/9). The resulting 
values of c 3 are presented in the second column (TPS/3) of Table |I[ We carried out the 



PA^ 


C 3 WAS*)'- 


PA/3 


c 3 (ECH): TPS/3 


PA/3 


c 3 (TPS-PMS): TPS/3 


PA/3 


[2/1] 


21.7 


21.7 


35.1 


35.1 


35.1 


35.1 


[3/1] 


13.7 


15.7 


19.5 


22.9 


19.0 


21.5 


[4/1] 


11.1 


15.8 


14.4 


20.8 


13.1 


18.7 


[5/1] 


9.3 


16.9 


11.2 


19.6 


8.8 17.3 


[1/2] 


12.8 


12.8 


17.3 


17.3 


17.3 


17.3 


[2/2] 


12.4 


14.9 


16.9 


20.4 


16.2 


19.4 


[3/2] 


11.7±3.4i 


15.8 


15.8±6.4i 


20.7±2.8i 


15.4±7.4i 


17.3±3.6i 


[4/2] 


10.3±2.8i 


15.7 


12.9±5.1i 


20.4±1.8i 


11.6±6.8i 


17.0±2.6i 


[1/3] 


12.4 


15.0 


16.9 


20.6 


16.2 


19.5 


[2/3] 


12.9 


15.1±1.2i 


17.4 


19.3 


18.3±0.8i 


18.5 


[3/3] 


10.6±2.9i 


14.0±1.7i 


13.6±5.5i 


20.2±2.0i 


12.6±7.0i 


16.9±2.7i 



TABLE I. Predictions for C3 in our, ECH and TPS-PMS approximants, using PA's of the invariant 
Borel transform B(z) of the approximants and demanding that the IRi pole be at z po i c = l//?o (= 4/9). 
"TPS/3" denotes that the parameters c£ ] {k > 4, J = 1,2) in and c k (k > 4) in ECH and TPS-PMS, 

are set equal to zero; "PA/3" denotes that the /3-functions in the approximants are resummed as: [2/3]^ 
(RSchl); [2/4]/3 (RSch2; cf ] =cf ) ); [3/2]^ (ECH RSch, and TPS-PMS RSch). 



analogous c 3 -fixing for the polynomial approximants ECH and TPS-PMS[] to the BjPSR, 
and c 3 predictions for them are also included in Table | (columns with "TPS/3"). These 
entries in the Table suggest the values c 3 ~ 12.5, 17, 16 for y/As?, ECH, and TPS-PMS, 
respectively. The predictions of PA^-'s of intermediate order ([3/1], [4/1], [2/2], [3/2], [1/3], 
[2/3]) appear to give the most stable results. Predictions of the higher order PA-g-'s gradually 
lose predictability (predicted c 3 's can even become complex) because of the afore-mentioned 
overdetermination. The lowest order PA^-'s are unreliable due to their too simple structure. 

The possibility to adjust the N 3 LO coefficient r 3 of (14]) in a similar way, was apparently 
first mentioned by the authors of Ref. f2"0fl . They referred to PA's ([2/1]) of the simple Borel 
transform, so their (PA-resummed) predictions would depend on the choice of the RScl 
and RSch. A systematic method to optimize the perturbative expansion by including the 
information on the location of the IRi pole was suggested in Ref. pl|. 



1 The ECH approximant is ^4^ CH ^(c 3 ) = a(lnQ| CH ; Pz? c 3; • • ■)'■> the TPS-PMS approximant is 

4 PMS '( c 3) = a PMS-P2apMs/ 2 ' with a PMs(c3) = a(lnQECH;3/52/2,C 3 ,...), QECH = ( 2o ex P(- r l/A))- 



Up until now we have taken the higher order parameters cjj, (k > 4, j — 1, 2) in our 
approximant (and in the ECH and TPS-PMS) to be zero, thus truncating the corresponding 
/3-functions (TPS/3). However, since the considered observable has low process energy Q p h ~ 
1.73 GeV, we expect the higher order terms occ^x k+2 (k > 4, x = a s /ir) of the /3-function 
to contribute significantly to the determination (via evolution) of the relevant coupling 
parameters of the approximants. This leads us immediately to the question of quasianalytic 
continuation of the (3(x) functions from the small-x into the large-x regime. We can choose 
again Pade approximants (PA's) as a tool of this quasianalytic continuation, keeping c 3 as 
the only free parameter, and subsequently determine c 3 in the afore-mentioned way. 

In y/As^ics) there are /3-functions characterized by the RSch-parameters (c 2 , c 3 , . . .) 

(2) 

(RSchl) and (c 2 , c 3 ,...) (RSch2) and determining the evolution and values of a x = 
a(\nQl; c 2 , c 3 , . . .) and a 2 = a(ln Q\\ c 2 , c 3 , . . .), respectively. In the (NNLO) ECH and 
the TPS-PMS approximants, the RSch-sets are (p 2 , c 3 , . . .) and (3p 2 /2, c 3 , . . .), respectively. 
For RSchl, ECH and TPS-PMS RSch, we have at first the freedom to construct [2/3], [3/2], 
or [4/1] PA^'s. For RSch2, the additional condition C4 = c 4 (5c 4 = 0) has to be fulfilled. 
Since C4 is a unique function of c 3 once a PA/31 choice has been made for RSchl, we then 
have for PA^ of RSch2 the possibilities [2/4], [3/3], [4/2], [5/1]. For each choice of PA^'s, we 
essentially repeat the afore-mentioned procedure of determining the value of c 3 . We consider 
the best choice of PA^'s the one giving the most stable prediction of c 3 over various PA's 
[M/N]-g of the approximant's invariant Borel transform. This turns out to be for V 'As 2 ( c 3) 
the choice ([2/3]^, [2/4]^), although ([2/3]^, [5/1]^) give virtually the same and almost 
as stable c 3 -predictions. For the ECH and TPS-PMS the choice is [3/2]^. The predictions 
for c 3 are given in Table | (columns with "PA/3"). Those from PA's [M/N]-g of intermediate 
order are significantly more stable than the corresponding ones with truncated /^-functions 
("TPS/3"). This is a numerical indication that the PA-resummation of the /3-functions im- 
proves the ability of the approximants to discern nonperturbative effects in the considered 
observable. The "PA/3" -entries in Table | give us approximate values c 3 = 15.5,20, 19 for 
our, the ECH, and the TPS-PMS approximant, respectively. 

There is yet another argument in favor of the above PA/3 choices. The chosen [2/3]/ji and 
[2/4]/32 (or: [5/1]^) have positive poles with mutually similar values: x po i e = 0.334,0.325 
(or: 0.291), respectively. The value of x po i e (=a vo i e /ii) indicates a point where "a strong and 



an asymptotically-free phase share a common infrared attractor" . Thus, it is reasonable 
to expect that only those RSch's whose /^(x)-functions have about the same value of x po \ e 
are suitable for the use in calculation of nonperturbative effects (on the other hand, in 
purely perturbative QCD, all RSch's are formally equivalent). Hence, the mutual proximity 
of Xp i e 's of RSchl and RSch2 PA/3's is now yet another indication that these PA/3's are the 
reasonable ones. What happens if we choose for RSchl and RSch2 other PA/3's? In such 
cases, we always end up with one of the following situations: Either the two corresponding 
positive Xpoie values are far apart, or both values are unphysically small, or one (positive) 
Xp ie doesn't exist, or there are no predictions for c 3 (not even unstable), or x po i e values 
are unstable under the change of c 3 in the interesting region c 3 ~ 12-16. So, the choice 
[2/3]/3i and [2/4] ^2 (or [5/1]^) in our approximant is not just the choice giving the most 
stable c 3 -predictions, it is also the only choice giving mutually similar (and reasonable) 
values of x po i e of RSchl and RSch2. Further, the choice [3/2] p for the ECH and TPS-PMS 
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RSch's gives us x po i e values similar to the ones previously mentioned: x po i e = 0.263 for ECH 
with C3 = 20; £ po i e = 0.327 for TPS-PMS with C3 = 19. Even other choices of PA^ for the 
ECH and TPS-PMS RSch's ([2/3]^, [4/1]^), which also give rather stable and very similar 
c 3 -predictions, give us x po i e ~ 0.27-0.41. Therefore, we see in all cases a clear correlation 
between the stability of the C3-predictions on the one hand and x po ie ~ 0.3-0.4 on the other 
hand. Finally, [2/3]^ is then the good choice for MS RSch since it has x po i e = 0.311 (c^ s , 
cf. Eq. (|i~5D, has been determined in Ref. [p3|| ). The choices [3/2] ^ and [4/1]^ for MS give 
Xp ie = 0.119, 0.213, respectively, which are further away from 0.3-0.4. 

Now that all the hitherto unknown parameters in \J~A of (H) and in the ECH and TPS- 
PMS have been determined, we use the approximants to predict the values of a^ s (3GeV 2 ) 
(=7ra ) from the measured values of the BjPSR 5 , (Q ph = 3GeV 2 ). Experimental values at 
Qph = V% GeV are given in [f23j] (their Table 4) and are based on SLAC data 

-\g A \ [l - S(Q 2 ph )] = 0.177 ±0.018 =>• S(Q 2 ph ) = 0.155 ± 0.086 . (19) 



where the constant \g A \ is known [24] from /3-decay measurements: \g A \ = 1.257 (±0.2%). 
The experimental uncertainties are high, mainly because of the effects of perturbative evo- 
lution on the small-XBj extrapolation of the polarized structure functions appearing in the 
sum rule (|IB|), as explained in Ref. We vary ao in our, and any other, approximant for 
the BjPSR S in such a way that the values (|19D are reproduced. We then obtain the predic- 
tions for a^ s (3GeV 2 ) given in Table O. Given are always three values for a s , corresponding 
to the three values of S (0). The results are given for our, the ECH and the TPS-PMS 
approximants, all with the described C3-fixing and with the afore-mentioned PA-type re- 
summation of the pertaining /3-functions: [2/3]p (RSchl), [2/4]^ (RSch2), [3/2]^ (ECH), 
[3/2]^ (TPS-PMS), [2/3]/3 (MS). Given are also predictions of such approximants when 
the /3-functions are TPS's (c^ = for k > 4). To highlight the importance of C3-fixing, 
we included predictions of these approximants (with TPS/?) when we set C3 = in them. 
In addition, predictions of the following approximants are included in Table 0: TPS Swi 
(0) (NNLO TPS); TPS S m with r 3 = 128.05 (N 3 LO TPS); off-diagonal Pade approximants 
(PA's) [1/2)5 and [2/l]s; square root of the diagonal PA (dPA) [2/2] 52, which is based solely 
on the TPS Sp] ([14]) , as are the previous two off-diagonal PA's; [2/2] 5 is the dPA constructed 
on the basis of the N 3 LO TPS S [3] with r 3 = 128.05. For [2/2] s and N 3 LO TPS we took the 
value r 3 = 128.05 (in MS, at RScl Qg = 3GeV 2 ) because then the [1/2] PA of the invariant 
Borel transform Bs (P"8j) predicts the correct IRi pole z po \ e = l/ j3 . Numbers in Table [H] are 
with four digits so that predictions of various methods can be easily compared. 

Table || includes predictions for a^ s (M|). They were obtained from af s (3GeV 2 ) by 
evolution via four-loop RGE, using the values of the four-loop MS coefficient c 3 (rif) f23]| 
and the corresponding three-loop matching conditions [^5| for the flavor thresholds. In the 
matching, we used the scale n(rif) = Km q (rif) above which rif flavors are active, with k = 2, 
and m q (nf) being the running quark mass m q (m q ) of the nj'th flavor. If increasing k from 
1.5 to 3, the predictions for a^ s (M|) decrease by at most 0.15%. If we use [2/3]^ instead 
of TPS/3 in the evolution from 3GeV 2 to Mf , ap(Mf) decreases by less than 0.04%. 

In Fig. |I] we present various approximants as functions of a^ s (M|). Our, ECH and 
TPS-PMS approximants have TPS /^-functions and c 3 = 12.5, 17, 16, respectively (by the 
described IRi pole requirement). These three approximants, when the /3-functions are 
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TABLE II. Predictions for a^ 18 , derived from various resummations of the BjPSR at Q^ h = 3GeV 2 . 
Predictions corresponding to S' max = 0.241 cannot be made with the TPS-PMS, because the latter cannot 
be larger than (2/3) 3 / 2 p 2 1 ' 2 fh 0.233, due to its specific polynomial form. 



resummed by PA's, and c 3 = 15.5,20,19, respectively (IRi pole), are presented in Fig. ||. 
The three approximants with TPS^'s (from Fig. [I]) are included in Fig. ^ for comparison. 

If we reexpand the approximants in powers of a (RScl Ql = Qpi l , in MS, n/ = 3), predic- 
tions for coefficient r 3 at Oq of expansion are obtained. Our approximant vl4.s 2 (c 3 ), with 
c 3 = 15.5, predicts r 3 = 125.8-cp/2+c 3 «130.8, and the ECH ^ ECH) = <z(lnQ| CH ; p 2 , c 3 , . . .), 
with c 3 = 20., predicts r 3 = 129.9 + (— c 3 vIS + c 3 )/2 « 129.4. Both predictions agree well with 
that of [^6j] r 3 « 129.9 (~ 130.) which was obtained from the ECH under the assumption 
(— c 3 4S +c 3 )~0 (note that c| IS ~21.0 was not known at the time [25] was written). 



IV. DISCUSSION AND CONCLUSIONS 



The results presented in Table [TI| and in Figs. [T]-0 show clearly that nonperturbative 
effects, as reflected in the mechanism of quasianalytic continuation from the small-a into 
large-a regime and in the presence of the leading infrared renormalon (IRi) pole, play an 
important role in the BjPSR at low photon transfer momenta Q p ^ ~ 1.73 GeV. These 
effects decrease the predicted value of a^ /[S (M|) by very substantial amounts. Our approx- 
imant gives the BjPSR-prediction ap(Mf) = 0.1120^;^ (see Table 0). Availability 
of additional data on polarized structure functions, especially in the Iow-xbj regime, may 
significantly reduce the uncertainties of the BjPSR-predictions for a^ s (M|). 

The present world average is a^(M|) = 0.1173±0.0020 by Ref. and 0.1184±0.0031 
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by Ref. pf. The NNLO TPS predictions of the considered BjPSR (O.imjftoajo, see Table 



n|) cover the entire world average interval and more. However, when the afore-mentioned two 
classes of nonperturbative effects are taken into account, e.g. via the use of our or the ECH 
approximants and by the described C3-fixing, we obtain an upper bound a^ s (M|) max « 
0.117 - see Table [HJ. But this upper bound does not surpass the central values of the 
afore-mentioned world averages. The central value of a^ s (M|) extracted from the BjPSR 
(f» 0.112) is significantly lower than the world average. 

What could be the reason for this? One speculative possibility would be that some of 
the Feynman diagrams contributing to the N 3 LO term (yet unknown) of the BjPSR have a 
genuinely new topology not appearing in the lower diagrams, and that such new topology 
diagrams push the predicted values of a^ s (M|) significantly upwards. The resummation 
methods based on the NNLO TPS cannot "foresee" such contributions IPjpBfl. In this context, 
we note that the afore-described C3-fixing in our, ECH and TPS-PMS approximants enables 
these approximants to be based on more than just the information contained in the NNLO 
TPS and in the RGE. However, since the location of the IRi pole can be determined by 
large-/3o considerations, the described C3-fixing apparently does not incorporate information 
from those possible higher-loop diagrams whose topologies are genuinely new. 

Another possible reason for the difference between our a^f -predictions and those of the 
world average could for example lie in a hitherto underestimated relevance of nonperturbative 
contributions and of higher order perturbative terms in the numerical analyses of data for 
some QCD observables. This possibility should be seen also in view of the fact that (some) 
NNLO contributions (~a 3 ) are not yet theoretically known for several of the quantities 
whose data have been analyzed to predict the world average [pTp8| . However, lower values 
are allowed by some recent analyses beyond the NLO: af s (M 2 z ) = 0.118 ± 0.006 |2§ from 
the CCFR data for x Bj F 3 structure function from uN DIS (NNLO); 0.112±g;gS| from Gross- 
Llewellyn-Smith sum rule |28| (NNLO); 0.115 ± 0.004 |27| from lattice computations. 



From the theoretical point of view, we are dealing with three types of resummation 
approximants for NNLO TPS's of QCD observables in the present paper: 

1. Pade approximants (PA's) provide an efficient mechanism of quasianalytic continua- 
tion. However, they do not possess RScl- and RSch-invariance, although their depen- 
dence on the RScl and on the leading RSch-parameter c 2 is in general weaker than 
that of the original TPS. In addition, they implicitly possess a c 3 -dependence, but this 
dependence has no special role since there is also c 2 - and RScl-dependence. 

2. Grunberg's ECH and Stevenson's TPS-PMS methods do not possess a strong mech- 
anism of quasianalytic continuation, except the one provided by the RGE-evolution 
of the coupling parameter a itself.Q This is so because these approximants do not 
go beyond the polynomial form in terms of the coupling parameter a. On the other 
hand, these approximants do achieve RScl- and C2-independence, since they represent 
a judicious choice of the RScl and of c 2 in the TPS. They possess a c 3 -dependence. 



In the one-loop limit, this amounts to the [1/1] PA quasianalytic continuation for a (ECH). 
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3. Our approximants provide an efficient mechanism of quasianalytic continuation, since 

1 /2 

they reduce to the diagonal PA expression [2/2] 02 in the one-loop limit (when all 
cjj. , cjj? , Cfc 1— ► for fc > 1). At the same time, they possess invariance under the change 
of the RScl and of the leading RSch-parameter C2. They possess a C3-dependence. 



The dependence on C3 (and on c^, k > 4) parameters in our, ECH and TPS-PMS 
approximants allows us to incorporate into them important nonperturbative infor- 
mation about the location of the leading IR renormalon pole. Further, it allows us 
to use in these approximants resummed /3-functions (PA-type), thus presumably ad- 
ditionally strengthening the effects of quasianalytic continuation mechanism. These 
approximants are then fully independent of the RScl and RSch of the original TPS. 



The leading higher-twist term contribution to the BjPSR (~1/Qp h ) P0|| - PT| , or a part 
of it, is implicitly contained in our approximant, as well as in the ECH and the TPS-PMS, 
via the afore-mentioned C3-fixing. The described approach implicitly gives an approximant- 
specific prescription for the elimination of the (leading IR) renormalon ambiguity. It is not 
clear which approximant accounts for the ~1/Qp h terms in the best way. 

A more detailed and extensive presentation of the subject will appear shortly f32|j . 
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FIG. 1. Values of various approximants (of Table [Tl|) as functions of (M§ ) ■ The /^-functions 
in our, ECH and TPS-PMS approximants have TPS form and the values of C3 in them were 
determined by the described IRi pole requirement. The experimental bounds ( |l9|) 5 m i n = 0.069, 
<S'max = 0.241 and <S m id = 0.155 are indicated as horizontal lines. 
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FIG. 2. As in Fig. [j], but the /3-functions in our, ECH and TPS-PMS approximants are now resummed 
by PA's and the values of C3 subsequently determined by the IRi pole requirement (see the text). For 
comparison, the corresponding predictions from Fig. [IJ with TPS/3 's, are included. 
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